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ABSTRACT 

We present a simple method for including the oblateness of a rapidly rotating neutron star when 
fitting X-ray light curves. In previous work we showed that the oblateness induced by rotation at 
frequencies above 300 Hz produces a geometric effect which needs to be accounted for when modelling 
light curves to extract constraints on the neutron star's mass and radius. In our model X-rays are 
emitted from the surface of an oblate neutron star and propagate to the observer along geodesies of 
the Schwarzschild metric for a spherical neutron star. Doppler effects due to rotation are added in the 
same manner as in the case of a spherical neutron star. We show that this model captures the most 
important effects due to the neutron star's rotation. We also explain how the geometric oblateness 
effect can rival the Doppler effect for some emission geometries. 

Subject headings: stars: neutron — stars: rotation — X-ray: binaries — relativity — pulsars: general 



1. INTRODUCTION 

The observation of pulsed X-ray emission originating 
from the surfaces of rotating neutron stars provides an 
excellent opportunity to study the properties of neutron 
stars and to possibly constrain the equation of state 
(EOS) of supernuclear density matter. Detailed mod- 
elling of the properties of the emitting region on the star 
combined with relativistic raytracing can be used to cre- 
ate theoretical light curves which can then be compared 
with the observed light curves to constrain the masses 
and radii (and hence the EOS) of the X-ray emitting 
neutron stars. 

The best candidates for a pulse-shape analysis are the 
accreting neutron stars where the energy released from 
accretion leads to X-ray emission from the surfaces of 
the neutron stars. The accreting neutron stars that are 
of interest fall into three classes: slowly rotating X-ray 
pulsars such as Her X-l ; accreting ms-Period neutron 
stars exhibiting Type I X-ray bursts; and the accret- 
ing ms-Period X-ray pulsars such as SAX J1808. 4-3658 . 
(These classes are not mutually exclusive: some of the 
accreting ms pulsars are also Type I bursters.) In the 
case o f Her X-l (with a spin period of 1.24 s), iLeahvl 
(|2004h has shown that an accretion-column model can 
explain the data and constrain the dimcnsionless mass- 
to-radius ratio to t h e ran ge 0.18 < M/R < 0.19. 
Bhattacha rvva et al.l ((2005) modelled the burst oscilla- 
tions detected during Type I X-ray bursts originating 
on XTE J1814-338 (spin period 3.2 ms) and found an 
upper limit of M/R < 0.24. An analysis of the X-ray 
pulse shape of SAX J1808. 4-3 658 (spin period 2.5 ms) by 
iPoutanen fc Gierlihskil (|2003h also provided constraints 
on the mass-to-radius ratios to the range 0.18 < M/R < 
0.3. In addition, one non-accreting ms X-ray pulsar 
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PSR J0437-4715 has been modelled by iBogdanov et all 
(2006) and has been shown to be consistent with a radius 
R > 6.7 km. While none of these constraints are strong 
enough to conclusively rule out any equations of state, 
these studies do show a great deal of promise. 

The most commonly used approximation used to 
model neutron star light curves is the Schwarzschild 
+ Doppler (S+D) approxim ation (|Miller fe Lamrj|1998t 
IPoutanen fe Gierlihskil l2003h . In the S+D approxima- 
tion the gravitational light-bending effects are mod- 
elled as though the star is not rotating using the 
Schwarzschild metric and th e form alism prescribed by 
iPechenick. Ftaclas. fe Cohenl (|1983h . Rotational effects 
are added by introducing Doppler terms as though the 
star is a rotating object with no gravitationa l field. In 
our previous calculations (jCadeau et a 1. 1 [2007: we found 
that as long as the star's spin frequency is below 300 
Hz, the ratio M/R can be accurately extracted using the 
S+D approximation. 

In two recent papers ((Cadeau et all 120051 l200l we 
investigated the validity of the S+D approximation by 
computing the geodesies of rapidly rotating neutron stars 
using the exact numerical metric to describe the gravi- 
tational field and th e shape of the star. We have shown 
(|Cadeau et al.ll2005h that it is necessary to correctly bin 
the light curve so that the variable times of arrival of 
photons from different parts of the star are taken into 
account if the star is rapidly rotating. However, we also 
showed that it was sufficient to use the Schwarzschild 
geometry to calculate the times of arrival: the effect of 
frame-dragging on light-bending and times of arrival is 
not large enough to affect t he light curves, co nsistent 
with the calculations do ne bylBraie et ail (120001 ). 

In the second paper ((Cadeau et al.ll2007l ) we investi- 
gated the effect of the neutron star's oblat e shape. Pre- 
vious calculations (Br aie fc Romanil |2001| ) of magneto- 
spheric scattering found that the oblate shape of a rotat- 
ing neutron star can alter the directions that photons are 
scattered into. In our computations (without scattering) 
we found that the oblate shape has a significant effect 
on the light curves which should be included when mod- 
elling the light curves of rapidly rotating neutron stars. 
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The reason for the importance of oblateness is simply 
geometrical: when light is emitted from an oblate sur- 
face the directions that the light can be emitted into are 
different from those possible from the surface of a spher- 
ical star. This leads to certain spot locations on the star 
where the spot is visible if the surface is oblate but would 
be invisible if the surface were spherical (and vice versa) . 

In our previous work (|Cadeau et al.1 120071 ) we found 
that as long as the correct shape of the star is used to 
formulate the initial conditions, the light curves result- 
ing from tracing rays in the Schwarzschild metric are a 
very good approximation to the light curves resulting 
from ray-tracing in the full numerical metric of a rapidly 
rotating neutron star. The problem with the method 
used in our past work is that it is too slow for use in fit- 
ting light curves to theoretical models. In our past work 
a full relativistic stellar structure calculation as well as 
geodesies connecting all spot and observer latitudes had 
to be computed for each hypothetical neutron star model. 
Since many thousands of models must be computed in or- 
der to fit the neutron star's properties to a light curve, 
our past work was not useful for comparing theoretical 
calculations with X-ray timing data. 

The purpose of this paper is to present a simple ap- 
proximation that includes the essential features of our 
more detailed ray-tracing calculations for rapidly rotat- 
ing neutron stars. Our new approximation, the Oblate 
Schwarzschild (OS) approximation, is based on the more 
commonly used S+D approximation. In the OS approx- 
imation an empirical formula for the oblate shape of 
the rotating neutron star is used to define initial con- 
ditions and visibility conditions for photons. Once ini- 
tial conditions on the surface of the star are set up, the 
Schwarzschild metric is used to find bending angles and 
photon time delays. Doppler effects due to rotation are 
then added to the model in the same manner as in the 
S+D approximation. 

The outline of this paper is as follows. In Section [5] 
we present a simple empirical formula for the shape of 
the star and derive expressions for the vector normal to 
the surface, angles between the surface normal and the 
light rays, and the solid angle subtended by an oblate 
surface. In Section [3] we explain how these concepts are 
used to create light curves in the OS approximation, pro- 
viding details of how the new visibility condition is im- 
plemented. In Section [4] we explain why the geometric 
effect of oblateness can rival the Doppler effect for some 
emission geometries. In Section [5] we compare our ap- 
proximate light curves with light curves computed using 
our previous exact methods. Finally, we conclude with 
comments on possible applications of the OS approxima- 
tion. 

2. THE OBLATE SHAPE OF A ROTATING NEUTRON STAR 

In the OS approximation, an oblate surface is embed- 
ded in the metric of a non- rotating neutron star. The 
metric exterior to a non-rotating neutron star is given 
by the Schwarzschild solution 

-l 



ds 2 



1 



dr 2 



+r 2 (d9 2 + sin 2 9d<p 2 ) . (1) 

(We use gravitational units where G = c = 1.) The 
surface of the star is described by a function r = R(9), 



and will be specified in subsection 12. H Given a function 
describing the surface of the star, the surface area of a 
small spot of angular extent dO and d<f) located at angles 
9 and <f> on the surface of the star will have a surface area 

dS(9) = R 2 {9) sin 0(l + / 2 (0)) 1/2 dBduft (2) 

where the function f(ff) is defined by 

(1-2M/R)- 1 / 2 dR (l + z) dR 
W) = E — = — (3) 



R 



d6 



where the gravitational redshift z = l/wl — 2M/R — 1 
has been introduced. 

An alternative (but equivalent) description of the grav- 
itational field exterior to a non-rotating star is given 
by the Isotropic Schwarzsch i ld me tric (see, for example, 
iMisner. Thorne. & Wheeler! (|1973f )) 



ds 2 = 



(1 - M/2f) 



di 2 



(1 + M/2f), 
+ (1 + M/2f) 4 (df 2 + f 2 (d9 2 + sin 2 9d<p 2 )) , (4) 

where the isotropic radial coordinate f is related to the 
areal radial coordinate r by the transformations 



and 



r = f (1 + M/2ff (5) 
dr = (1 - 2M/r) 1/2 (1 + M/2f) 2 df. (6) 



The surface of the star in this coordinate system is de- 
noted f = R(6). Making use of the coordinate transfor- 
mations given by equations ([5|) and ([6]) it is easily shown 
that the function f(9) defined in equation ([3]) is 



1 dR 

R~dJ' 



(7) 



2.1. A Model for Oblateness 



In order to find the shape of a rotating neutron star, 
it is necessary to specify an equation of state (EOS) of 
dense nuclear matter, a spin period and mass and then 
solve the relativistic equations of stellar structure. The 
relativistic stellar structure equations must be solved 
numerically using an axisymm etric code such as rns 4 
(|Stergioulas fc Friedman! I1995D. based on the m ethods 
described by iCook. Shapiro, fc Teukolskvl ((l99l . Once 
the stellar structure equations have been solved, the lo- 
cation of the star's surface is extracted by finding the 
locations where the fluid's enthalpy vanishes. 

In order to aid the process of fitting light curves, we 
have found a simple empirical formula for the shapes of 
rotating neutron stars. We assume that the surface of a 
rotating neutron star can be described by a function of 
the form 



R{6) 



R, 



= 1 



i-q 



2 

n=0 



(C,e)P2„(cos0) 



(8) 



where 9 is the co-latitude, measured from the spin axis 
(cos (9 = on the equator), P n (cos9) is the Legendre 
polynomial of order n, R eq is the radius of the rotating 

4 Available at http://www.gravity.phys.uwm.edu/rns 
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star measured at the equator and the parameters £ and 
e are given by 

C = ^ (9) 



and 



GM 



c 2 C 



(10) 



with f2 = 2ir/P, where P is the spin period. The coef- 
ficients a n (C,e) can then be fit to the sh apes of stellar 
models computed using the rns code. The lHartle] (|1967t ) 
slow rotation approximation corresponds to dropping the 
<Z4 term in the series. We find that the 04 term is required 
to describe the shape of the largest, most rapidly rotating 
stars. However, higher order terms are not required. 

We have computed the stellar structure of compact 
stars for a wide variety of equations of state. For each 
EOS, models spanning the allowed values of mass and 
angular velocity were computed. We find that the EOS 
can be grouped into two families, where each family's 
stars can be described by a separate set of Legendre co- 
efficients. The first family, Neutron a nd Hybrid Quark 
Stars, includes the following EOS: lArnett fc Bowera 
(11973) catalogue EOS A , B, C, F, G, L, N, O; EOS 
APR (|Akmal et alj I1998T ) which includes nuclei scat- 
tering data and special relativistic corrections; hybrid 
qu ark stars with a mixe d quark-hadron phase calculated 
by Alford et al.l (120051); and Hype ron stars computed 
by lLackev. Navvar. fc Owenl (20061 ) using methods de- 
scribed by iGlendennind (|2000 ). The second family con- 



sists of color-flavor-locked (CFL) quark stars described 
by the EOS P = (e-4B)/3 where B is the bag constant. 
Bag constants in the range 50 — 140 MeV/fm 3 are used 
to generate di fferent CFL EOS. This simple EOS has 
been shown bv lFraga et al.l (|2001h to be an excellent ap- 
proximation to quark star models including second order 
pcrturbative QCD corrections. 

Although each EOS has a different mass versus radius 
curves, for any given combination of mass, radius and 
angular velocity the scaled shape of the star is indepen- 
dent of the EOS within its family. The values of the 
Legendre fitting coefficients a n ((,e) for the two families 
are given in Table 1. Since R(tv/2) is defined to be the 
equatorial radius of the rotating star, this requires that 
the coefficients obey = a — a 2 /2 + 3a 4 /8. As can be 
seen from the coefficients in Table 1 this requirement is 
not exactly satisfied. However, the errors in this model 
are still less than 1%. 

Equation © is only useful if the equatorial radius of 
the star is known. When doing light curve fits, we will 
need to choose a value of 9 describing the spot's co- 
latitude and the value of radius at the spot's latitude. 
Given these values, equation ||5J) can be inverted to solve 
for iL„, 



R, 



< </ 



R{6) 



= l + ^ 6 2™(C,e)^2„(cos( 



n=0 



+P 2 (cos#)^c 2n (e~)P 2 „(cos#), (11) 

n=l 

where C = GM/R(9)c 2 and e = Q 2 R 3 (9)/GM. The val- 
ues of the b n and c n coefficients are displayed in Table 1. 



It is our intention that equation (jlip be used in light 
curve fitting to put constraints on the equatorial radius 
of a rotating star. In order to compare with the predic- 
tions of a specific EOS, it is still necessary to construct a 
sequence of neutron stars with a specific spin frequency. 
The construction of such a sequence for a fixed spin fre- 
quency can be done in a straight-forward fashion using 
a public domain program such as rns. 

2.2. The Normal Vector 

The isotropic form of the Schwarzschild metric given in 
equation iJH has the useful feature that surfaces of con- 
stant time are conformally flat. This allows the intro- 
duction of quasi-Cartesian coordinates {x,y 7 z} defined 
by x = fsm9cos(f>, y — fsin#sin0 and z = fcos9. 
In the quasi-Cartesian coordinate system the z axis is 
aligned with the star's spin axis, and the observer lies in 
the x-z plane at an angle i from the spin axis. 

The vector normal to the surface at any arbitrary point 
can be constructed using the quasi-Cartesian coordinate 
system in the same way that it would be computed in flat 
space. First a tangent vector is constructed by taking the 
gradient of the function R making use of equation ([7]). 
The normal to the surface at the location of the spot is 

n = (1 + f{9)) ~ 1/2 [(sin 9 - cos 9f(9)) (cos <j) x + sin y) 

+ (cos0 + sin 0/(0)) z] . (12) 

The unit radial vector is defined by r = sin 9 (cos x + 
sinewy) + cos 9 z, so the angle 7 between r and n, is 
simply 



cos 



7=(l + / 2 (0)) 



-1/2 



(13) 



Since sin 7 = //yl + / 2 , the normal vector in equation 
(fT2|) can also be expressed as 

n = sin(0 — 7) (cos (j> x + sin <f> y) + cos(0 — 7) z (14) 

as would be expected from simple geometric considera- 
tions in flat space. 

2.3. Bending Angles and Zenith Angles 

A light ray emitted in an initial direction 1 will be mov- 
ing a different direction k when it reaches the observer. 
Since the observer is located at an angle i from the spin 
axis, k = sini x + cosi z. In the standard treatment of 
spherical neutron stars, a is the angle between the ini- 
tial photon direction and the radial direction, defined by 
cos a = 1 • r and the bending angle tjj is defined by 

cos ip = k • r = cos i cos 9 + sin i sin 9 cos cf>. (15) 

In the Schwarzschild metric, the vectors 1 and k are 
coplanar so that 

1 = . - , (sin ak + sm(ib — a)r) . (16) 
sm ip 

The zenith angle (3 between the normal vector and the 
initial photon direction is defined by cos/3 = 1-n. Making 
use of equation (fTB)) . the identity n = (sin(0 — 7)r + 
sin 7 z)/ sin 9, and the spherical trigonometric identity 
cos i — cos 9 cos -0 + sin 9 sin if; cos 8 (see Figure [1] for the 
geometrical definition of 5) the zenith angle has the value 

cos/3 = 1 • n = cos a cos 7 + sin a sin 7 cos S, (17) 

if sin tp ^ 0. If sinV> = 0, it is straight forward to show 
that cos P = cos a cos 7. 
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Fig. 1. — Spherical triangle illustrating the emission geometry. 
The vector z is aligned with the spin axis, r points in the direction 
of the spot, and k points in the direction of the observer. The 
angles between these three vectors as well as the interior angles of 
the spherical triangle are shown. 

2.4. Solid Angle for Static Oblate Stars 

In this section we derive the solid angle subtended by 
an surface area element on an oblate star that is at rest. 
In the subsequent section we will discuss the modifica- 
tions necessary when the star is rotating. 

Since the photons travel to the observer via geodesies 
of the Schwarzschild metric, the solid angle subtended 
by a surface element is given by the standard expression 
(|Pechenick. Ftaclas. & Cohenlll983l ) 



bdb. 



D 2 



(18) 



where b is the photon's impact parameter, D is the dis- 
tance between the star and the observer, and (p is the 
azimuthal angle around the vector k. The subscript "s" 
refers to static, since in this section we do not take into 
account Doppler shifts and boosts due to the rotation of 
the star. The impact parameter is related to the zenith 
angle a by 

b f 2M~ 
s [ na = ——Jl - (19) 



R(9) 



R(ey 



The bending angle tp is related to the impact parameter 
and the surface of the star by 



iP(b,R(8)) 



R(9) 



dr 

~2 



1 



r(l 



2M 
r 



-1/2 



(20) 



In order to express the solid angle in terms of the co- 
ordinates defined on the surface of the star, a change 
of variables from coordinates {6, 4> } to {?/>, <j>o\ must be 
made so that 



db 



bdipd(j) 
LP ' 



(21) 



The relationship between the impact parameter and the 
bending angle for a spherical star is simple, since the pho- 
tons are emitted from a surface of constant gravitational 
potential. In the case of an oblate star, the bending angle 



ip depends on both the impact parameter and the vari- 
able location of the star's surface. The partial derivative 
appearing in equation (|2ip is 



(db_ 



db\ fdb\ 
di>) R + \dR) y 



dR 

dip 



(22) 



where {db/dip) R = 1/ {dip/db) R is evaluated using the 
integral definition given in equation (|20|) . Making use of 
standard relations for partial derivatives 



dR), 



db\ fdjA 
chl>) R \dR) b 



(23) 



Using equations (fT9f and (|20|) the second partial deriva- 
tive appearing on the right-hand side of equation (123| 
is 

fchjA (1 + 2:) sina 

\dRj b R cos a' [ ' 

The partial derivative of R with respect to ip can be 
calculated using equations §5§ and (flT 



1 (^L 

R Uhp 



f(0) 



(1 + z) sinff \ dip 



dcosO 



(25) 



Making use of standard identities for the spherical trian- 
gle defined in Figure [TJ 



d cos 

dip 



sin 8 cos S. 



(26) 



The cos S term in this last expression can be eliminated 
using equation (|17[) . leading to the simple equation 



R \~chp 



(cos (3 — cos a cos 7) 
(1 + z) cos 7 sin a 



The final result for equation 



db\ (db\ 



(27) 



(28) 



The expression for the solid angle can now be simplified 



to 



dfl s - (l + z ) 2 ^cos(3 



d cos a 



dcosip 



simp 
cos 7 



dipdepo. (29) 



The transformation between the {ip, <po\ coordinates to 
the star's coordinates {9, 0} is given by the usual relation 
sin ipdipd(p = sm9d0d(p. Equation ([2]) can now be used 
to simplify the solid angle formula to 



dS 

dQ s = (1 + z) 2 — ^ cos/3 
D z 



d cos a 



dcosip 



(30) 



The solid angle subtended by a spot on an oblate star 
given by equation (|30[) is very similar to the expression 
for the solid angle subtended by a spot on a spherical 
surface. The expressions for the oblate and spherical sur- 
faces differ through the cos (3 term, corresponding to the 
cosine of the angle between the initial photon direction 
and the vector normal to the surface. In the spherical 
case, cos/3 is replaced by cos a, the cosine of the angle 
between the initial photon direction and the radial vec- 
tor. Hence our result for the solid angle reduces to a 
simple correction for the tilt of the surface induced by 
oblateness. 
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3. LIGHT CURVES FOR ROTATING OBLATE STARS 

The flux due to a spot on the star is 

dF = I dfl (31) 

where I is specific intensity measured in the observer's 
frame and d£l is the solid angle subtended by the spot 
as measured by the observer. 

Light emitted with frequency v em will detected with 
frequency v a given by 

V 



1 + z 

where 77 is the Doppler boost factor given by 



V = 



vT 



(32) 



(33) 



1 — v cos £ ' 

where v is the magnitude of the spot's velocity and £ is 
the angle between the velocity and the original photon 
direction. The spot's velocity is unaffected by the oblate 
shape of the star, so the functions v and £ are given by 
the values 

-1/2 



v = 2Trv*R(0)sin6 1 - 



2M 
~R 



and 



cos£ = — - 



sm a sin 1 sin < 



sin-0 



(34) 



(35) 



The motion of the spot causes an aberration in the 
angles between the normal to the surface and the initial 
photon direction as measured in the rest frame of the 
spot, j3 em and the observer at infinity, /3, given by 

cos (3 em — 77 cos (3 (36) 

where (3 is given by equation (jTTJ) . In order for the spot 
to be visible, it is necessary that < cos/3 em < 1. The 
solid angle measured at infinity dil is also affected by the 
spot's motion, and is related to the solid angle derived 
in the previous section by 

dtt a = ndfl s . (37) 

A photon emitted from R{6) and with impact param- 
eter b will arrive at a time 

T(b,R(6)) = f°° dr (l-— } x 
Jr(9) \ r J 

relative to a photon with zero impact parameter emitted 
from the same location. 

Finally, making using of the relativistic conservation 
law I /^o = I l v emi the measured flux is 



F„ 



dn s r] 4 (l + z) d I(casf3 em ,i> em ) (39) 



COS /3em>0 



where dQ. s is given by equation (|30[) . In order to pro- 
duce a light curve the flux should be integrated over the 
observed frequency range, and the flux should be binned 
according to the relative arrival times of the photons. 
Times of arrival can be computed exactly using equa- 
tion (I38p or a series expansion such as that given by 
iPoutanen &: Beloborodovi (|2006l ) (see their equation 18) 
can be used. 
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Fig. 2. — Side view of a rotating neutron star (solid curve). The 
surface of a spherical star with the same radius as the oblate star 
at the emission latitude is shown as a dashed circle. A photon with 
initial direction 1 is initially moving towards smaller values of the 
radial coordinate. At r c the photon reaches the minimum value 
of r and begins to move outwards. The change in bending angle 
while the photon travels between the two intersections of the solid 
and dashed curves is 2Aip. 

3.1. Initially Ingoing Photons 

For photons that are emitted in an outgoing direction, 
the calculation of the flux is straightforward and similar 
to the standard calculation for a star with a spherical 
surface. However, because the surface is oblate, there will 
be photons emitted which are initially directed radially 
inwards. These photons will travel to smaller values of r 
until a critical radial coordinate r c is reached, at which 
point r = 0. After the critical radius has been reached, 
the photon will move outwards. In Figure [5] an initially 
ingoing photon's trajectory is labelled by vector 1. 

An initially ingoing photon will have a value of impact 
parameter 

b in < . R = , (40) 
y/l - 2M/R 

and an angle a > ir/2. The quantity R is the location of 
the surface of the star at the latitude of the spot. The 
critical radius is given by the solution of the equation 



1 



2M 



(41) 



Given this critical radius, the bending angle associated 
with the motion of the photon from R to r c is Aif) given 
by the integral 



rf?" 



2M 



-1/2 



(42) 



Similarly, the relative coordinate time taken by the pho- 
ton to travel this distance is 

r R , / 2A/V 1 ( Hn ( 2M^^~ 1 ' 2 



(43) 

By symmetry, when the photon travels from r c out to R 
the bending angle is Alp and the extra time is AT. The 
final result is that the total bending angle for an initially 
ingoing photon is 



^in (h n ,R) = 2 Alp + 1p(bi n , R) 



(44) 
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where ip(b,R) is given by equation ((20)) . Similarly, the 
relative arrival time for the photon is 



T in (b in ,R) = 2AT + T(b in ,R). 



(45) 



As expected, an initially ingoing photon will have a larger 
bending angle and take longer to reach the observer than 
an initially outgoing photon with the same value of im- 
pact parameter. 

In all situations of interest, the difference R{9) — r c 
will always be much smaller than the radius of the star. 
This makes it possible to approximate the integrals Atp 
and AT by assuming that in the integrand r = r c (l + e) 
where e <C 1. To lowest order in e 1 / 2 , 



A*p = 



, (R-r c ) 
' (r c - 3M) ' 



(46) 



(47) 



Similarly an approximation for AT is 

AT = r c (l - 2M/r c )- 1/2 Aif;. 

4. EFFECT OF OBLATENESS ON LIGHT CURVES 

It may seem counter-intuitive that the oblatcness of a 
rotating star could have an important effect on the re- 
sulting light curves since a naive argument predicts that 
oblateness is unimportant: since oblateness is an effect 
that is second-order in angular velocity it ought to be 
smaller than hrst-order effects such as Doppler effects 
and t he frame-dragging effect. However, ICadeau et al.1 
(2007) have shown that the effect of oblateness can be 
large in some cases. The main reason why oblateness 
can't be neglected is not through its effect on the star's 
gravitational field. It is instead a simpler geometric argu- 
ment that would also hold in flat space. We now explain 
why the naive argument fails for certain ranges of spot 
locations. 

Doppler shifts and boosts are first-order corrections to 
the shape of a light curve, so Doppler corrections are the 
most im portant effect. As has been sho wn elsewhere (for 
example iPoutanen fc Gierlm ski (2003)) Doppler boost- 
ing creates an asymmetry in the light curve. However, 
Doppler effects are most important when the spot is mov- 
ing in the same direction as the initial photon direction. 
For light emitted near the back of the star (for values 
of rotational phase near (p — tt) the angle £ defined in 
equation (|35[) is close to tt/2 so that the boost factor 
in equation (|33[) only differs from unity by corrections 
that are second order in v/c. It is in this region near 
the back of the star that oblateness effects become most 
important. 

Frame-dragging is almost always a more significant 
contribution to the gravitational field than the contri- 
butions to the field due to oblateness. Only at very high 
spin frequencies when e ^ 0.1 (corresponding to about 
500 Hz, depending on the EOS) do the oblateness cor- 
rections to the spacetime metric become comparable to 
the frame-dragging corrections. However, neither correc- 
tion has a strong influence on the trajectories of photons 
or the resulting light curves. As the photon moves away 
from the star, the frame-dragging frequency falls off as 
r -3 so it rapidly becomes unimportant. (Corrections due 
to oblat eness fall off even fa ster.) In our previous calcu- 
lations (jCadeau et al.l [20071 ) we found very little differ- 
ence between the light curves computed with the Kerr or 



Scwarzschild metric. For this reason, the use of the Kerr 
metric (which keeps the frame-dragging term in the met- 
ric) is not a great improvement over the Schwarzschild 
metric (which ignores frame-dragging) when computing 
light curves. 

Oblateness affects lights curves in two ways: the solid 
angle is changed since the emitting surface is tilted at 
a different angle with respect to the initial photon di- 
rection; and the visibility condition is altered. The first 
effect is usually small (except near the back of the star) 
while the second effect can be surprisingly large. 

The first oblateness effect is easily seen in equation 
(f30|) for the solid angle of an oblate star. This expression 
is equivalent to the expression for the solid angle of a 
spherical star multiplied by the factor cos /3/ cos a. This 
factor corrects for the change in the tilt of the surface 
due to oblateness. Since cos/3 (the angle between the 
surface normal and the initial photon) only differs from 
cos a (the angle between the radial direction and the ini- 
tial photon) by an amount of order 0(f2 2 ) this correction 
is not typically large. However, near the back of the star 
(in relation to the observer) this correction factor can 
be important if cos a is of similar order of magnitude 
as 7 or smaller. In order to see this, consider the value 
of phase <p> = tt, corresponding to the back of the star. 
At this point, equation (fT7|) for the zenith angle reduces 
to cos (3 — cos a cos 7 + sin a sin 7. If cos a is small and 
O(cosa) ~ 0(7) then cos/3/ cosa ~ 1+7/ cosa4-0(7 2 ). 
In this case the correction factor cos (3 / cos a is not nec- 
essarily close to unity, and can cause large changes in the 
received flux. Similar arguments hold for other values of 
phase as long as cos a is small. 

The second effect of oblateness, due to the change in 
the visibility condition, can be understood with a sim- 
ple example. Consider a spot at fixed latitude 9 on two 
equal-mass stars, a spherical star and an oblate star, cho- 
sen so that the radii of the two stars at the spots' latitude 
are the same. Suppose that at the value of rotational 
phase denoted (pi the spot is at the limb of the spherical 
star. This means that a = tt/2 and the bending an- 
gle is at its maximum value, ipmax- Using equation (I15|) 
cos ipmax = cos i cos 9 + sin i sin 9 cos (pi . 

Similarly, on the oblate star, the spot is at the limb at 
rotational phase <j>2 = 4>\ + A(p at which point (3 = tt/2 
and the angle between the radial direction and the initial 
photon direction is ct2, differing by tt/2 by an amount sec- 
ond order in angular velocity. The bending angle for light 
emitted from <f>2 is cos ip2 = cos i cos 9 + sin i sin 9 cos (f>2 ■ 
The difference between the rotational phases on the two 
stars where the spot is at the limb is given by 

COS 1p2 — COS Ipmax ( d COS 1p \ COS «2 

COS <p2 — COS (pi = 



sin i sin f 



dcosa /„ sin i sm( 



( 48 ) 

If we assume that the phase difference, A(p is very small, 
and that the spot is very close to the back of the star so 
that <pi ~ tt, then cos <p2 — cos (pi — (A(p) 2 /2. As a result, 
the change of phase over which the spot is visible due to 
oblateness is approximately 

/ \ 1/2 

A0~ — — - V -\ cosa 2 49 

\ sinismfl \dcosa/ Q2 J 

when the spot is near the back of the star. In the 
iBeloborodovl ((20021 approximation d cos ip/d cos a ~ (1+ 
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z) 2 for light emitted close to the normal to the surface. 
In our example light is emitted close to the tangent to the 
surface, but the cosine derivative continues to be finite 
and order unity. Since cos 012 ~ 0(VL 2 ) it follows that 
the phase difference A<fi ~ 0(Q). As a result, for light 
emitted on the side of the star opposite the observer, the 
phase differences caused by oblateness are first order in 
the angular velocity. Since the Doppler boosting is only 
second order in the angular velocity in this region of the 
star, oblateness becomes the dominant effect. 

5. COMPARISON OF MODELS 

The effect of oblateness is most important for very 
large neutron stars rotating with very high spin frequen- 
cies. In order to test the accuracy of the OS approxima- 
tion discussed in this paper, we now examine the results 
of this approximation for one very oblate neutron star 
with a mass of 1.4M Q , an equatorial radius of R eq = 16.4 
km and a spin frequency of 600 Hz. The expansion pa- 
rameters defined in equations ([9]) and (fT0|) take on values 
of C = 0.126 and e = 0.339 for this model. 

Consider first the case of emission from a spot at co- 
latitude 9 = 49°. The radius of the star (computed using 
the exact numerical method) at this point is R(8) = 15.1 
km. The oblateness model given by equation JHJ) dif- 
fers from this value only in the fourth significant figure. 
In Figure [3] light curves for light emitted from this co- 
latitude and observed at an inclination of i = 70° are 
shown. The bold solid curve labelled "exact" is com- 
puted using the ex act numerical method described by 
ICadeau et all (|2007ft . The curve labelled "S+D" is com- 
puted for a spherical 1.4M Q star with a radius of 15.1 
km. The curve labelled "OS" is computed using the 
Oblate Schwarzschild approximation described in this 
paper. For this combination of emission and observer an- 
gles, the S+D approximation predicts an eclipse for part 
of the rotational period. In the cases of the exact calcu- 
lation and the OS approximation the spot is visible all of 
the time. The small differences between the Exact and 
OS light curves are due to differences in the spacetime 
metric and to small errors due to the oblateness model 
given in equation (JSj> . The differences between the OS 
and Exact light curves are very small compared to the 
differences between the S+D and Exact light curves. 

In Figure 2] light curves for the same star are shown, 
but in this case light is emitted from a co-latitude of 
8 = 41°, and the observer's inclination is i = 20°. For 
this stellar model, the radius of the star at 8 — 41° is 
14.8 km. The light curves in Figure [4] have the same 
meaning as in the previous figure, however for the S+D 
approximation, the surface of the spherical star is cho- 
sen to be 14.8 km, in order to keep the values of the 
gravitational potentials the same in both OS and S+D 
approximations. As in the case for emission from 49°, the 
OS and Exact light curves differ by an amount which is 
much smaller than the difference between the S+D and 
the Exact light curves. 

6. CONCLUSIONS 

In our previous work (|Cadeau et al.l [20071 ) we demon- 
strated the need to include the oblateness of rapidly ro- 
tating neutron star when computing light curves. In 
this paper we presented a simple method, the Oblate 
Schwarzschild approximation which can be used to ac- 




FlG. 3. — Light curves for light emitted from a co-latitude of 
9 = 49° and an observer at an inclination of i = 70° . The neutron 
star has M = 1.4Mq, R eq = 16.4 km and spins at 600 Hz. The 
solid bold curv e is computed using t he exact numerical space time, 
as described bv lCadeau et all (120071 ). The solid curve labelled S+D 
is computed by assuming that the star's surface is spherical. The 
dashed curve is computed using the Oblate Schwarzschild approx- 
imation described in this paper. 
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Fig. 4. — Light curves for light emitted from a co-latitude of 
6 = 41° and an observer at an inclination of i = 20°. The neutron 
star's parameters are exactly as given in Figure [3] 

count for the star's oblate shape. The method pre- 
sented in sections [2] and [3] can be summarized as follows. 
First the mass, radius and angular velocity of the neu- 
tron star are chosen. The shape of the star's surface 
is then given by the empirical formula given in equa- 
tion ((SJ), which in turn allows the calculation of the an- 
gles 7 and /?, where 7 is the angle between the surface 
normal and the radial direction, and (3 is the angle be- 
tween the initial photon direction and the surface nor- 
mal. Given the angle j3, the light curve for the rapidly 
rotating neutron star can be computed using the same 
formulae us ed in the S+D approximation (as given, for 
example by iPoutanen fc Gierlihskil (|2003l )) with the fol- 
lowing changes: the angle a (between the initial photon 
direction and the radial direction) is replaced by /3; the 
visibility condition is given by < cos j3 < 1; and the 
photons must be binned by their times of arrival. In the 
OS approximation the bending angles and times of ar- 
rival are calculated using the Schwarzschild metric. The 
resulting light curves closely approximate the exact light 
curves, as shown in section [5] 

The accreting ms X-ray pulsars such as 

SAX J1808.4-3658 (|Wiinands fc van der Klisl Il99l 



TABLE 1 

Neutron Star Shape Parameters 



Coefficient 


Neutron and Hybrid Quark Stars 


CFL Quark Stars 


U() 


-0.18e + 0.23Ce-0.05e 2 


-0.26e + 0.50Ce - 0.04e 2 


U-2 


-0.39e + 0.29Ce + 0.13e 2 


-0.53e + 0.85Ce + 0.06e 2 


04 


+0.04e-0.15Ce + 0.07e 2 


+0.02e-0.14Ce + 0.09e 2 


bo 


+0.18?-0.23C? + 0.18f 2 


+0.26e-0.50Ce + 0.31? 2 


In 


+0.39e-0.29Ce + 0.42? 2 


+0.53e-0.85^e+ 1.06? 2 


hi 


-0.04f + 0.15C?-0.13? 2 


-0.02e + 0.14£?-0.12? 2 


C2 


+0.60? 2 


+1.13e 2 


C4, 


-0.12? 2 


-0.07? 2 



IChakrabartv &; Morganl |1998[ ) are rotating rapidly 
enough that light curve models should include the 
oblateness model presented in this paper along with 
realistic models of the emission spectrum and geometry. 
Comparisons of the light curves arising from these oblate 
models w ith data for SAX J 1808. 4-3658 are currently in 
progress (|Leahv et al.ll2007f) . 
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